On the two-steps relaxation of mean-field glasses: p-spin model 
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Critical slowing down dynamics of supercooled glass-forming liquids is usually understood at the 
mean-field level in the framework of Mode Coupling Theory, providing a two-time relajcation scenario 
and power-law behaviors of the time correlation function at dynamic criticality. In this work we 
derive critical slowing down exponents of spin-glass models undergoing discontinuous transitions 
by computing their Gibbs free energy and connecting the dynamic behavior to static "in-state" 
properties. Both the spherical and Ising versions are considered and, in the simpler spherical case, a 
generalization to arbitrary schematic Mode Coupling kernels is presented. Comparison with dynamic 
results available in literature is performed. Analytical predictions for the Ising case are provided for 
any p. 



I. INTRODUCTION 

The slowing down of the dynamics of supercooled 
glass-forming hquids corresponds to a non trivial under- 
lying thermodynamic landscape. An unusual time be- 
havior of the density correlation function emerges, with 
respect to the exponential decay: a separation between 
fast (/3) and slow (a) relaxation modes takes place and 
the correlation function develops a plateau approaching 
a dynamic arrest transition. ^^'^ 

Theoretical advances have been made studying such 
systems with the so called Mode Coupling Theory 
(MCT),^"® a mean- field theoretical description of many 
particle systems, able to identify the separation of two 
relaxation processes and a dynamical transition tempera- 
ture Td at which ergodicity breaks down with the system 
undergoing a structural arrest. In this framework, for 
T ^ Td the time spent at the plateau increases and di- 
verges at Td- At criticality a power-law behavior arises 
for the correlation close to the plateau and the exponents 
governing the approach to t~°') and the depart from 
(~ t^) the plateau are related as: 



r(l-2a) ^ r(l + 2fe) 



(1) 



where A is a model dependent quantity (functional of the 
static structure factor^), that is usually treated like a 
tunable parameter. 

Since the works of Kirkpatrick, Thirumalai and 
Wolynes,^""^^ the behavior of glass-forming liquids and 
structural glasses has been linked with dynamics and 
thermodynamics of a certain class of mean-field spin- 
glass (SG) models, sometimes called mean-field glasses 
or discontinuous spin-glasses. These include the p-spin 
models, ^^^^^ displaying a dynamic transition at which 
ergodicity breaks down, with a behavior of the corre- 
lation function close or identical to the one predicted 
by MCT. Below this dynamic transition the Boltzmann 
measure splits in many well defined metastable states, 
which at a temperature become thermodynamically 



stable inducing a phase transition. ^-'^'^^'-'^^ The static tran- 
sition is governed by an entropy crisis and plays the 
role of the hypothesized Kauzmann transition, related 
to the vanishing of the configurational entropy of the 
liquid. Below this point, named Random First Or- 
der Transition (RFOT), the model develops a spin-glass 
phase with one step replica symmetry breaking (IRSB). 
Besides the conjectured thermodynamic analogy between 
discontinuous spin-glasses and structural glasses, cf. Mo- 
saic Theory, ^^'^^'^^ the former can as well be exploited 
for studying critical slowing down dynamics and MCT 
properties at the dynamic arrest transition. 

In this paper we do not address to the study of the 
thermodynamic but we investigate only the dynamical 
transition: we apply to the p-spin model a recent devel- 
opment in the comprehension of the dynamical transition 
and its relation to the thermodynamic behavior, i.e., a 
relationship between critical slowing down exponents and 
thermodynamic "in-state" quantities:^''^''^^ 



Wi 



(2) 



with the A defined in Eq. (1). 

The coefficients wi, W2 are obtained from the expan- 
sion at the third order terms of the appropriate Gibbs 
free-energy r[g] functional of the order parameters, i.e. 
the replica overlap matrix elements Qap . The potential F 
is the Legendre transform the free energy functional as 
function of replica "pinning" fields. It coincides with 
the Franz-Parisi potential, by means of which the dy- 
namic transition is identified as the spinodal point of an 
excited local minimum at a non-zero overlap value equal 
to the plateau value of the correlation function. Critical 
slowing down exponents can, then, be obtained from the 
coefficients of the expansion of the F potential around a 
replica symmetric (RS) solution q'^^ at the plateau value 
of the correlation function at the dynamic critical point. 

For T ^ Td the Gibbs potential can be expanded in 
powers of overlap fluctuations around the dynamic solu- 
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tion, Sqafj = Qaf} - Q^f: 

6T[Sqal3] ~ ^ Mai3,^S SqafjSq^S + 
(a/3), (75) 

Wai3f5e<p ^qafi^fl'i&^fle'f , (3) 
(a/3)(75)(£¥') 

where ST[Sqa0] = T[Qap] - ^Qap] 2ind the first order 
term is absent at the saddle-point. Here and below brack- 
ets in sums over replica indices mean that only distinct 
replicas must be considered. 

The dynamic transition is associated to the vanishing 
of one eigenvalue of the mass term, the replicon, which 
defines a critical direction in the replica space. Projecting 
this expansion along the replicon direction S^^^qab (see 
next session for the details) one obtains: 

SnS^'-haff] ~ «;iTr(5Mg)3 + W2J2 (^^'^ea/?)" . (4) 

In order to obtain the MCT exponents for mean-field 
discontinuous SG models one can proceed with the fol- 
lowing protocol: 

1. compute the averaged replicated action T[Q]; 

2. compute the expansion around the dynamical RS 
solution up to the third order, cf. Eq. (4); 

3. solve the system equation for the saddle point 
and the vanishing of the replicon, both in the 
n — > 1 limit which allows to work on the dynamic 
metastable state of the Gibbs potential; 

4. evaluate the coefficients wi, W2 of the third order 
along the replicon direction; 

5. compute A by means of Eq. (2) and a and h with 
Eq. (1). 



II. THE SPHERICAL P-SPIN MODEL 

In this section we focus on the spherical version of the 
fully-connected p-spin model. Its relevance in the con- 
text of mean-field glasses is due to the fact that for tem- 
peratures down to the dynamic transition, the dynamic 
equations for its two time correlation functions^^'^^ are 
equivalent to those of the schematic MCT.^'^ The Hamil- 
tonian of the model reads: 

H = - ^ ...(Tip - /i^CTi , (5) 



where the couplings J are Gaussian independent identi- 
cally distributed variables with: 



P(J) 



7rp! 



■ exp 



-J' 



(6) 



and the spins a are real values variables subject to a 
global constraint: 



(7) 



Due to the spherical constraint this model is analytically 
solvable in all details needed for our scope. This allow us 
to show how to compute analytically third order coeffi- 
cients and dynamic exponents a and b step by step. 

Through a saddle point calculation, it is possible to 
compute the replicated partition function:^^ 



Zn 



„S(cx>) 



/ X{\l^dq^pe^v{-NV[Q]] 
J\q\>o „-^w-, V 27r 



'l9l>0a</J 

■2pZ^^«/3 2 



(8) 



where the overbar means average over disorder, cf. Eq. 
(6), M = Qa0 = <^i<^i and 



^^(oo) 



gAr(l-Mog(27r))/2^-l/2 



l + O 



N 



(9) 



Differentiating Eq. (8) with respect to Qa0 one gets the 
saddle point condition: 

mQ^^' + + (O-')a^ = , (10) 

where the term proportional to /i^ is absent because it is 
irrilevant in both the n — ^ 0, 1 limits. 

For our purpose we also need the third order Taylor 
expansion of F with respect to fluctuations around the 
RS solution: Sq^p = Qafi - O^f, with = (1 - q)d^ + 
g,where is the Kronecker delta: 



(a/3)(75) 



+ 4 E ^a0'ySev:>^<lcl3SqjsSqsv , (H) 



3! 

where Sr[Sq] = T[Q] - F[g^'5]. The first order is absent 
at the saddle point and 



-{p - i)/xQ:/<52;<5^ + {Q-'U{Q-')50 + m 



(12) 

4 



- a/37(5e(/? 



OQafjdQysdQeip 



(13) 
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In order to study the critical dynamic behavior we work 
within a RS Ansatz with n — >■ 126-28, 30 g^^^^^ restrict our 
analysis to the replicon subspace, defined by the condi- 
tions 











(14) 



(the second condition is a consequence of the first one as 
far is a symmetric matrix. The vanishing of the 

eigenvalue of the Hessian in this subspace, the so-called 
replicon eigenvalue, yields the criticality condition. 

We now restrict the analysis to the case without ex- 
ternal magnetic field h = Q. Imposing the saddle point 
condition, cf. Eq. 10, at the dynamic transition point 
(RS with n ^ 1) 



q 



(1-9) 

and the vanishing of the replicon 

1 



(p-1)M^-2 



(1-9) 



2 ' 



(15) 



(16) 



leads to the expressions for the value of the overlap at 
the dynamic transition and the dynamic temperature: 



Qd = 



p-2 
p-1 



(p-2)f-2 



(17) 



The dynamic transition point is also the point at which 
the dynamic saddle point solution appears, that is, Eqs. 
(15)-(16) are not independent and, indeed, Eq. (16) can 
be obtained also as the derivative of the Eq. (15). 

Considering fluctuation only in the replicon subspace, 
one can considerably simplify the expansion of Eq. (11) 
as 



25r[<5gW] -(p-l)(p-2)M^-3^(5g('-j)3 



{l-qf 



(18) 



(a/37) 



In this case the tensorial form of Eq. (13) is so simple 

that we can straightforwardly compute the values of the 
cumulants Wi and yielding 



W2 



{p-l){p-2) 



M^-3(l - qf 



which, imposing Eq. (16), reduces to: 

ip-m-i) 



W2 



2q 



(19) 



(20) 



Using the value of qd at the transition, cf Eq. (17), one 
obtains for A the p-independent value 



1 

2 ' 



(21) 



which coincides with the result reported in Ref. [18,31]. 



A. The case of uniform magnetic field 

The presence of a magnetic field term {h ^ 0) can 
change the nature of the transition. As shown in Rcfs. 
[17,18] for values of the field h> htr = ^p^-v{p- 2)^12 
the transition becomes continues: no plateau is there 
and at the transition the long time limit of the correla- 
tion function does not jump discontinuously, though the 
relaxation behavior in time is still a power-law in the ^ 
regime. The exponent a (sometimes called v) is the only 
one defined. We flrst study this continuous transition 
and then wo move to the discontinuous one. 

The value of a can be computed as for the discontinu- 
ous transition case, except for the fact that now one has 
to work in the n ^ limit. ^'^'^'^ This implies that the RS 
expression of Eq. (10) becomes 



(l-# 



(22) 



rather than Eq. (15). Equations (16) and (19), instead, 
do not change. From Ref. [17] we know that the transi- 
tion line is parametrically defined as 



j,2 _ 



p{p - 1) 



(i-9)V-' 



2 ^ p(p-2) 1 
2 ^ 



(23) 

(24) 



with 1 — 2/p<q'<l. Given a value of the field h, one 
can, thus, straightforwardly compute the corresponding 
values of q and T from Eq. (23) and (24) and obtain A 
from Eq. (20). 

For the discontinuous transition in a fleld h < htr, two 
non-zero overlap values arc relevant: the plateau value 
qi and the long-time limit qo = C{oo). Comparing with 
the results of Ref. [17] we observe that the exponent 
parameter A is still given by the ratio W2/w\, cf. Eq. 
(20) , now evaluated on g = gi , solution of the following 
equations for the dynamic critical values of Td(/i), qo{h) 
and qi{h): 



= l,{p-l)qr\ cf. Eq. (16) 



(1-ai)' 



and 



1-qi 1-qo 

p-i 



■ p—l p 
M ( <7i -% 



') 



which are the IRSB saddle point equations (cf. Eq. (10)) 



Generalization to arbitrary schematic MCT 
models 



Comparing the dynamics of the p-spin spherical 

model^^ with the MCT differential equation for the corre- 
lation function, ^'^ one can notice that the function iJLq^~^ 
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{ficf)P~^ , in the MCT notation), derivative of the first term 
of the action (8), plays the role of the MCT memory 
kernel. One can generalize this argument, through 
schematic MCT, to a generic polynomial kernel^ 

Aim = n{v}, m = J2 vp^Pitr-'- (25) 



We start considering the long-time limit of the correlation 
function, (else called non-ergodicity parameter /), and 
the (p{t) expansion around such limit: 

ct>{t)c^qd + G{t); « 1 , (26) 

whose Laplace transform reads 

<^(z) ~ ~^+G{z) . (27) 



These behaviors have to satisfy the MCT dynamical 
equation and its transformed: 

Tdt^{t) + <j){t) = - I duA{t - u)du(j>{u) (28) 

<A(z) 



1 + Z(j){z) 



IT + K{z) . 



(29) 



Plugging Eq. (27) into the left hand side of Eq. (29) and 
expanding up to the third order, one obtains 



-qd+zG{z) 
l-qd+zG(z) 



(j>{z) 
^l+z4>{z) ^ 

l-S + fg#+0(G3)l . (30) 



The Laplace transform of the expansion of the memory 
kernel reads: 



zA{z) ~ -F{v,qd) + F'{v,qd)G{z) 

+ lF'\v,q,)CT[Gm{z) 



(31) 



One has to equate equations (30) and (31) order by order 
in G{t). The zeroth and first order are the standard long- 
time MCT equation and its derivative, cf. Eqs. (15), (16): 



F{v,qd) = 
F'{v,qd) = 



1 - 

1 



(l-?c/)^ ■ 

The third order yields: 

zG{zf + \zCT[G{t)\z) = Q 

with: 

A= V>,5d)(l-gd)3 . 



(32) 
(33) 

(34) 
(35) 



Assuming a power-law solution G{t) ^ {to/t)"' ^ one gets 
back Eq. (1) with a parameter exponent coinciding with 
Eq. (35). 



This result can also be obtained studying a model 
whose action is a slight generalization of Eq. (8) action: 



r[Q] 



2 A' 



A(Q„^)-ilog|0| (36) 



such that A(x) = A'{x). For polynomial function A{x) 
this action describe a model with an Hamiltonian com- 
posed by a proper sum of p-spin interaction terms, like 
Eq. (5).^^^^* Prom the action Eq. (36) one can eas- 
ily derive Eqs. (32), (33), with a kernel given by Eq. 
(25), yielding critical temperature (identical to the mode 
coupling temperature) and critical plateau value of the 
correlation. Expanding Eq. (36) to third order, cf. Eq. 
(11), yields: 



A=^ = iA"(,,)(l-g,)3 



(37) 



coinciding with Eq. (35) and verifying the method pro- 
posed for the schematic MCT models. This is the sim- 
plest case, where equations of motion can be solved and 
analytic results are available. We now move to consider 
a more difficult case for which A cannot be computed 
directly solving the dynamic equations and Eq. (2) re- 
mains the only way, known so far, to estimate the critical 
slowing down exponents. 



III. THE ISING P-SPIN MODEL 

In this section we focus on the Ising version of the fully- 
connected p-spin model. The Hamiltonian of the model 
reads: 



n = - 



■CTi. 



(38) 



where the couplings J are Gaussian distributed again 
with Eq. (6) and cxi = ±1. It is well known^^ that this 
model displays first a RFOT with a dynamical transition 
Td separated from a static transition to a spin-glass 
IRSB stable phase and, at a lower temperature, a sec- 
ond transition to a spin-glass full RSB phase. For our 
purpose, we only focus on the first transition. 

Averaging over disorder and introducing the overlap 
matrix through an auxiliary matrix A^b, one gets the 
replicated action: 



r[Q,A] = 

(a/3) 



2 ^ 

(a/3) 



AapQafi 



+ logTr{<,} W[A;ct] 



>V[A;£7] = exp j i ^ AapaaCrp | 
V («/3) / 



(39) 



(40) 



that has to be evaluated through a saddle point calcula- 
tion. The derivative with respect to A^^ yields 



(41) 
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leading to 

r[Q] 



(«/3) 



W[Q;a] = exp ( ^ ^ Qa/jVacr/j 

(a/3) 



(42) 



This concludes the first step of our protocol. The first 
order derivative reads: 



r' — 



dT 



dQ 



(< aa<Ji3 > -Qafi) , (43) 



where 6 = p{p — l)/?^/4 and < • • • > means average over 
the weight W[Q], of. Eq. (42). The vanishing of this 
equation gives the saddle point condition, that, solved in 
the RS Ansatz with n — ^ 1, yields 



< >= Af-'^ j dz W[z] tanh(^)2 
tanh(i;) 

exp ( - ..a9^..^-i I cosh(2:) 



q 

m 

W\z] — ^..j^ , , 

/p — 1 
dz W[z] = ■s/irpPq~ exp 

The second order derivative reads: 



(44) 

(45) 



(46) 



where the presence of connected averages (. • .)c is 
a consequence of the direct derivative of the term 
logTr^r W[Q;fT] (cf. App. A). In order to impose crit- 
icality, F"^ evaluated at the saddle point condition 
and projected onto the replicon subspace, should vanish. 
The first part of Eq. (46) is proportional to F^^, cf. Eq. 
(43) and it does not contribute. The vanishing of the 
second part, in the RS Ansatz, reads 



20qP-^{l -2q + r)~l = 



(47) 



as detailed in App. A. Here r =< m!^ >. This allows to 
rewrite Eq. (47) as: 



= OqP-^ < (1 - m2)2 >= 
= eqP-'^ < sech^(a;) > . 



(48) 



Once the saddle point and the vanishing of the second 
order derivative (Eq. (46)) are imposed, the third order 
derivative reads: 



/06-3p 



-pi/i 

03 al3,jS,E<p 



(p - 2)Q3'2p , 

\y J^afi X7 X<5 X V 
02 "a"fi'^a"l3 

+ < CFoi(J0(T^(TsCF£<J^p >c (49) 



p 


Ta 


Id 


A 


a 


b 


2 


1 





1/2 


.395 


1 


2.05 


.916 


.051 


.556 


.379 


.892 


2.2 


.808 


.198 


.652 


.346 


.768 


2.5 


.724 


.428 


.719 


.320 


.609 


3 


.682 


.643 


.743 


.308 


.570 


4 


.678 


.815 


.746 


.307 


.565 


5 


.700 


.881 


.743 


.308 


.570 


6 


.727 


.915 


.739 


.310 


.576 


1 


.loo 


.aoo 


. / OO 


.oil 


.051 


8 


.784 


.948 


.733 


.313 


.586 


9 


.812 


.957 


.731 


.314 


.589 


— >■ oo 


/ " 


1 


2/3 


.340 


.700 



TABLE I: Dynamic exponents in the Ising p-spin model. 



and this allows to write the coefficients of the expansion, 
cf. Eq. (4), as 



8613 



6-3p 
Id 

89^ 



Wi 



W2 = 



l-3qd + ivd -Ud = 
< (1 - m^f > 

^{qd - 2rd + Ud) + A = 

2 < 771^(1 - TTl^f > +A 



2{p - 2Y^-^P 



pyip-iy 



(50) 

(51) 
(52) 



where u =< >. As it happens in the Sherrington- 
Kirkpatrick (SK) model^^'"^^, the term A vanishes if p = 2 
and, indeed, can be considered as the correction to W2 
due to the multi-body interaction. 

With this results one can obtain the mimcrical values 
of the coefficients and the exponent a reported in Tab. L 



A. The p^2 limit 

The interest to the behavior of the model for p close 
to two is due to its relation with the SK model. Do- 
ing an expansion for small e = p — 2, one can expect 
and, actually, self-consistently verify that the finite jump 
of the overlap at the transition is of order e. One can 
consequently expand the action for small Qa/S, still con- 
sidering the transition as discontinuous. 
Prom Eq. (41) one obtains 



1 



^l7 = Qc 



(53) 
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Putting this result in the action (39) and expanding for 
small h-o^B one gets: 



+ iTrA3 + 0(A4) . 



2 

a/3 



a/3 



(54) 



The first three derivatives read 
2 



2r' 



a/3 



9r" 



or'" 



-A 



l-e 



2/32 
(2-g)(l -g) 

2^2 



+ (A2)„^ (55) 



asm 



(2 -£)(!-£)£ 



2/32 



(56) 



(57) 



From the first two derivatives one gets the criticality con- 
dition. In the RS Ansatz and with n — >■ 1 this reduces to 
the system equation: 



■A2 

2A 



where A is the off-diagonal part of A^ 
equation is solved by: 



Arf 
/3d 



f — ^elog£ . 



(58) 
(59) 

The system 

(60) 
(61) 



Evaluating the third order derivative (57) on this result 
one gets: 



A = 



W2 



(62) 



This result agrees with the a = .395 proposed by Kirk- 
patrick and Thinmialai^^ studying the dynamics of a 
soft-spin vcirsion of the model in the p — >■ 2 limit. 

We note that the behaviour of the parameter A is dis- 
continuous as a function of p at p = 2. Indeed in the SK 
model, i.e. precisely at p = 2, we have X — W2 — 0, while 
as soon as e > we have A = 1/2. This happens because 
the coefficient IV2 is proportional to the third derivative 
of (j2+'s which is singular at g = as soon as e is different 
from zero. 



B. The p — > 00 limit 

The behavior of the p-spin model for large p has been 
previously studied, due to its relation to the Random 



Energy Model (REM).'*'' In this section we, indeed, 
present the calculation of A = ■u;2/'*i'i in the limit p — 00. 

As in the p ^ 2 case, here is more convenient to work 
with the auxiliary variable A, the RS off-diagonal element 
of the matrix A^jS, which is related to q by Eq. (41). 
Furthermore in order to keep finite q, one should expect 
and consistently verifies that A diverges in the large p 
limit. From Eq. (44) one gets: 



a(A) 




=< l-(l-m2) > 



~ l-e-^/2\/2^, A»l 



(63) 



where, in the last expression, only the leading term, for 
A » 1, has been retained. Differentiating this equation, 
one arrives to the system equation for the criticality con- 
dition. At the leading order it reads: 



/(A) = 1-6-^/2^2^ 



2g(A) 
pX 



(64) 
(65) 



Solving this system first for A and then for q, which con- 
tains the whole dependence from the temperature, one 
finds: 



Ad ~ 21ogp 
Td ^ 



P 



41ogp 



(66) 
(67) 



As p — )■ 00 the critical dynamic temperature diverges and 
the dynamic overlap g'd — )• 1, cf. Eq. (64). 

Now we move to the computation of the exponent pa- 
rameter. Eq. (50), (51) and their ratio A, in the large A 
limit, reduce to: 



Wi 



W2 




A = 



Evaluating the last term one obtains: 



(68) 

(69) 

(70) 



Ae^/24 



/32A 



1 



277 log p 



1 



p y 41ogp 



-2p 



(71) 



which goes to zero as p — )• 00. 

Our result is, indeed, A — )• 2/3, corresponding to an 
exponent a = .340 
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C. The addition of ferromagnetic couplings 

It is possible to generalize the previous results allowing 
the couplings J to have a non zero mean: 



P{J) 



INP- 



Trpi 



■ exp 



NP- 



, (72) 



To treat the case Jo 7^ one has to introduce a non-zero 
magnetization for the system and the the Gibbs effective 
action, cf. Eq (42), becomes: 

T[Q, A,m,x] = ^J2 Qab + 13 Jo E <+ 

ab a 



ab 



+ logTr{,} W[A,x;a] 

W[A, X; ct] = exp f i E Kbf^af^b + ^ Xa(Ja J , 
\ ab a J 

(73) 

where the fields Xa play for the magnetization, rrio, the 
same role of A for the overlap. Through a saddle point 
calculation one arrives to two coupled equations: 

q = (m2)=A/'~^ / W[2]tanh(z)2 (74) 



m 



X 



(m) =J\f-^ j dz W[z] tanh(z) 
{z-ppJomP-^'}'^^ 



exp - 



/ 



PP'^Q' 



dz Wlz] = y/npPq~ 



(75) 



(76) 
(77) 



The solution with m = 0, the SG phase, is always present 

but, for Jo large enough, a ferromagnetic (FM) RS solu- 
tion appears discontinuously with m = q > and turns 
out to be the stable one. 

Since in the action (73) Jq couples only to the magne- 
tization m, the SG phase, where m = 0, is not afi^ected by 
the presence of a non-zero mean of the couplings. As a re- 
sult, along the whole dynamic SG transition the physics 
does not change and the values of the exponents a and b 
are constant. 

The PM/FM transition is a usual thermodynamic first 
order transition with a ferromagnetic spinodal line.^^ 

Two relevant point arc, indeed, present: the tricritical 
point between the SG, the FM and the FM phases and 
the intersection between the dynamic transition line and 
the FM spinodal line. It has been shown'^^''"^ that both 
the relevant points belong to the Nishimori Line (NL) 



(78) 



a line in the Jo , T phase diagram. This fact is due to the 
following property of systems on the NL:^° 



lim r(/3, 4 

n— >-0 



NL 



{P),q,m) 



limr(/3,0,q,0) 

n— >1 



(79) 



Equation (79) corresponds to the statement that the 
static Gibbs free energy potential along the NL is equal 
to the dynamic Gibbs free energy along the Jo = axis. 
This means that, in order to obtain the ratio W2/W1 at 
the dynamical transition with n — >■ 1, one can, simplify- 
ing the numerical calculation, evaluate those coefficients 
along the NL, at the spinodal point and working with 
n — >■ 0. Furthermore, since on the NL the melting pro- 
cess is equivalent to a glassy transition, one can argue 
that the exponent calculated through our procedure con- 
trols this melting process too. 



IV. CONCLUSIONS 

We have applied a new method^* to compute the slow- 
ing down exponents on the Ising and spherical versions 
of the frustrated ;)-spin model with Gaussian interaction. 
This method allows us to derive the mode coupling expo- 
nents of the critical dynamics by means of an analytical 
static-driven computation. Those exponents govern the 
power-law approach to {C{t) ~ + t~°') and departure 
from {C{t) ~ — t'') the plateau value qd. 

For the spherical case we exactly reproduce the analyt- 
ical result obtained in Ref. [18], derived by a Langevin 
description of the dynamics and, in full generality, all 
schematic MCT exponents. 

For the Ising version we present our computation and 
the exact values of both exponents for any value of p. 
In this case, the dynamics with discrete spins has not 
been directly solved so far and ours arc the first analytic 
estimates for a and b. One can, anyway, approximate 
discrete by soft spins in order to construct a dynamical 
equation. In that case the final computation for p > 
2 differs: the discrete case depends on the value of p, 
whereas the soft one does not. Our results agree with 
the computation performed in Ref. [12] on a soft-spin 
approximation in the p — )■ 2 limit. In this limit we found 
a discontinuity with the Sherrington-Kirkpatrick model 
{p = 2), with a finite jump, from to 1/2 of the exponent 
parameter A. For the sake of completeness, also the p 
00 limit has been characterized. The study of the Ising 
p-spin model concludes considering a non-zero mean of 
the couplings (Jq 7^ 0): up to a critical value, where a 
ferromagnetic transition takes place, the behavior of the 
model does not change with Jq, validating our estimates 
for the exponents along the whole SG transition line. 
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Appendix A: Computation of the cumulants in the 
replicon sub-space 

In this appendix we present the calculation of the con- 
nected cumulants of four and six replicas projected in the 
replicon sub-space. At first should be noticed that in the 
RS Ansatz four replica index quantities could take only 
three different values, depending on how many replica 
indeces are repeated. For: 

^(a/3)(75) =< <^a(Tl3CT^as >c = 

=< (Joi<y pcr^crs > — < OoiO^ >< u^a^ > , (Al) 
we have: 

= 1-9^ (A2) 

=9-9' (A3) 
^fa)(..) =r-<l\ (A4) 

where here a, /3, 7, 5 are considered different. This allows 
us to express the tensorial form of C^'^^ : 

+(r - e) . (A5) 

Summing over all replica indeces the C^^^ times the fluc- 
tuations in the replicon sub-space ( Eq. (14) ), the only 
term that does not vanish is (1 — 2g -|- r), with a factor 2 



due to the exchange of a with /3: 

^(a/3)(75)°^a/3"57'5 ~ 
(a/3), (75) 

= 2(1 -2g + r)^ (5^3)2 (A6) 

(a/3) 

Six replica index quantities, like: 

<^(a/3)(75)(e¥>) = < OaafiO-^a 1,(7^(7^ >c = (A7) 

= < aacr^a^asasa^ > 

— < aaCrpajas >< o^Gip > 

— < (J-yCFsCFe'^ip X CaC/3 > 

— < cFa(7 |3C^sa^p >< criers > 

+ 2 < (TaO'p >< cr^crg >< a^atp > , 

can take eight different values. The computation follows 
as in the previous case^^ leading to: 

(a/37) 

+16(g-2r + «)^(^gS)^ (AS) 

(a/3) 

which concludes the calculation. 
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